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\ 1 INTRODUCTION 
• i— I . 

■ Observations of temperature fluctuations in the cosmic microwave 
5^ ' background (CMB) provide a valuable means of distinguishing be- 
tween two competing theories for the formation of structure in the 
early Universe. Currently, the most favoured theory is the simple 
inflationary cold-dark-matter (CDM) model, for which the distri- 
bution of temperature fluctuations in the CMB should be Gaussian. 
The second class of theories invokes the formation of topologi- 
cal defects such as cosmic strings, monopoles or textures, which 
should imprint some non-Gaussian features in the CMB (Bouchet, 
Bennett & Stebbins 1988; Turok 1996). Thus, the detection (or oth- 
erwise) of a non-Gaussian signal in the CMB is an important means 
of discriminating between these two classes of theory. 

In order to test for large-scale non-Gaussianity in the CMB, 
the 4-year COBE-DMR dataset (in various forms) has already been 
analysed using a number of different statistical techniques, as dis- 
cussed below. These tests have been performed either on some 
combination of the 31-, 53- and 90-GHz A & B 4-year DMR maps, 
or the 4-year DMR maps from which Galactic emission has been 
removed. Two such Galaxy-removed maps are generally available, 
each one created using a different separation method. The DMR- 
DCMB map is a linear combination of all six individual COBE- 
DMR maps designed to cancel the free-free emission (Bennett et 
al. 1992), whereas the DMR-DSMB map is constructed by first 
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We investigate the detection of non-Gaussianity in the 4-year COBE data reported by Pando, 
Valls-Gabaud & Fang (1998), using a technique based on the discrete wavelet transform. 
Their analysis was performed on the two DMR faces centred on the North and South Galactic 
poles respectively, using the Daubechies 4 wavelet basis. We show that these results depend 
critically on the orientation of the data, and so should be treated with caution. For two distinct 
orientations of the data, we calculate unbiased estimates of the skewness, kurtosis and scale- 
scale correlation of the corresponding wavelet coefficients in all of the available scale domains 
of the transform. We obtain several detections of non-Gaussianity in the DMR-DSMB map at 
greater than the 99 per cent confidence level, but most of these occur on pixel-pixel scales and 
are therefore not cosmological in origin. Indeed, after removing all multipoles beyond I — 
40 from the COBE maps, only one robust detection remains. Moreover, using Monte-Carlo 
simulations, we find that the probability of obtaining such a detection by chance is 0.59. We 
repeat the analysis for the 53+90 GHz coadded COBE map. In this case, after removing I > 40 
multipoles, two non-Gaussian detections at the 99 per cent level remain. Nevertheless, again 
using Monte-Carlo simulations, we find that the probability of obtaining two such detections 
by chance is 0.28. Thus, we conclude the wavelet technique does not yield strong evidence 
for non-Gaussianity of cosmological origin in the 4-year COBE data. 
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subtracting templates of synchrotron and dust emission and then 
removing free-free emission (Bennett et al. 1994). 

The first investigation of non-Gaussianity in the 4-year COBE 
data was performed by Kogut et al. (1996). This analysis used the 
4-year DMR 53 GHz (A +B)/2 map at high latitudes (\b\ > 20°) 
with cut-outs near Ophiuchus and Orion (Bennett et al. 1996), and 
found that traditional statistics such as the three-point correlation 
function, the genus and the extrema correlation function, were com- 
pletely consistent with a Gaussian CMB signal. Colley, Gott & Park 
(1996) also computed the genus statistic, but for the DMR-DCMB 
map with \b\ > 30°, and arrived at similar conclusions. The full set 
of Minkowski functionals were computed for the 4-year 53 GHz 
(A + B)/2 map (with a smoothed Galactic cut) by Schmalzing & 
Gorski (1998), taking proper account of the curvature of the celes- 
tial sphere. They also concluded that the CMB is consistent with 
a Gaussian random field on degree scales. On computing the bi- 
spectrum of the 4-year COBE data, Heavens (1998) also found no 
evidence for non-Gaussianity. Finally, Novikov, Feldman & Shan- 
darin (1999) have calculated the partial Minkowski functionals for 
both the DMR-DCMB and DMR-DSMB maps and do report de- 
tections of non-Gaussianity, but the analysis was performed with- 
out making a Galactic cut and the detections most probably result 
from residual Galactic contamination. 

Recently, however, two apparently robust detections of non- 
Gaussianity in the 4-year COBE data have been reported. Ferreira, 
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Magueijo & Gorski (1998) applied a technique based on the nor- 
malised bi-spectrum to a map created by averaging the 53A, 53B, 
90A and 90B 4-year COBE-DMR channels (each weighted ac- 
cording to the inverse of its noise variance) and then applying the 
extended Galactic cut of Banday et al. (1997) and Bennett et al. 
(1996). They concluded that Gaussianity can be rejected at the 98 
per cent confidence level, with the dominant non-Gaussian signal 
concentrated near the multipole i = 16. This non-Gaussian signal is 
certainly present in the COBE data, but Banday, Zaroubi & Gorski 
(1999) have now shown that it is not cosmological in origin and 
is most likely the result of an observational artefact. Nevertheless, 
using an extended bi-spectrum analysis, Magueijo (1999) reports 
a new non-Gaussian signal above the 97 per cent level, even after 
removing the observational artefacts discovered by Banday et al. 

A second detection of non-Gaussianity was reported by 
Pando, Valls-Gabaud & Fang (1998) (hereinafter PVF), who ap- 
plied a technique based on the discrete wavelet transform (DWT) 
to Face and Face 5 of the QuadCube pixelisation of the DMR- 
DCMB and DMR-DSMB maps in Galactic coordinates (i.e. the 
North and South Galactic pole regions respectively). PVF com- 
puted the skewness, kurtosis and scale-scale correlation of the 
wavelet coefficients of DMR maps in certain domains of the 
wavelet transform, and compared these statistics with the corre- 
sponding probability distributions computed from 1000 realisations 
of simulated COBE observations of a Gaussian CMB sky. In all 
cases, they found that the skewness and kurtosis of the wavelet 
coefficients were consistent with a Gaussian CMB signal. On the 
other hand, the scale-scale correlation coefficients showed evidence 
for non-Gaussianity at the 99 per cent confidence level on scales 
of 1 1-22 degrees in Face of both the DMR-DCMB and DMR- 
DSMB maps. Nevertheless, in both maps, Face 5 was found to 
be consistent with Gaussianity. We note that Bromley & Tegmark 
(1999) confirm the findings of both PVF and Ferreira et al. (1998). 

In this paper, we also apply to the 4-year COBE data a non- 
Gaussianity test based on the skewness, kurtosis and scale-scale 
correlation of the wavelet coefficients. In the analysis presented 
below, however, we calculate the skewness and kurtosis statistics 
using unbiased estimators based on ^-statistics (Hobson, Jones & 
Lasenby 1999 - hereinafter HJL), as opposed to the straightforward 
calculation of sample moments employed by PVF. For the scale- 
scale correlation, we adopt the same definition as that used by PVF. 
We also note that the analysis presented below is slightly more gen- 
eral than that presented by PVF, since we calculate the statistics of 
the wavelet coefficients in all the available domains of the wavelet 
transform, as opposed to using only those regions that represent 
structure in the maps on the same scale in the horizontal and verti- 
cal directions. 

Perhaps the most important point addressed in the analysis 
presented here, however, is the fact that non-Gaussianity tests based 
on any orthogonal compactly-supported wavelet decomposition are 
sensitive to the orientation of the input map. This is discussed in de- 
tail below. As an illustration of this point, we therefore present the 
results of two separate analyses, in which the relative orientations 
of the input maps differ by 180 degrees. Nevertheless, it should be 
remembered that, in general, different techniques for detecting non- 
Gaussianity are each sensitive to different ways in which the data 
may be non-Gaussian. We should therefore not be too surprised if 
the detailed results of an analysis are orientation dependent. Obvi- 
ously, it would be troubling if the general conclusions concerning 
non-Gaussianity of the data depended on orientation, but that is not 
the case here. 



2 THE WAVELET DECOMPOSITION 

The basics of the wavelet non-Gaussianity test are discussed in 
detail in HJL and also by PVF and so we give only a brief out- 
line here. The two-dimensional discrete wavelet transform (DWT) 
(Daubechies 1992, Press et al. 1994) performs the decomposition 
of a planar digitised image of size 2 Jl x 2 Jl into the sum of a set of 
two-dimensional planar (digitised) wavelet basis functions 

7i-iy 2 -12 j i-12'2-l 

yW = E E E E b juhMuh^h,h4uhi x i)- (d 
1 ji=oj 2 =o /i=o ; 2 =o 

In equation (|l|), the wavelets Wji,j2\h,k( x ) ( w i m jlij2ih>h tak- 
ing the values indicated in the summations) form a complete and 
orthogonal set of basis functions. Each two-dimensional wavelet 
is simply the direct tensor product of the corresponding one- 
dimensional wavelets A|fi,;/, (■*) and *(/j 2; / 2 (y), which in turn are de- 
fined in terms of the dilations and translations of some mother 
wavelet via 

/2M 1/2 

J V^x/L-h), (2) 

where ^ x ^ L, and a similar expression holds for v/j 2 -j 2 (>')• Thus, 
the scale indices j\ and 72 correspond to the scales L/2 J1 anAL/2 11 
in the x- and y-directions respectively (so J\ and J2 are the small- 
est possible scales - i.e. one pixel - in each direction), whereas 
the location indices l\ and I2 correspond to the (x,y) -position 
(Ll\/2 ]l ,Ll2/2 J1 ) in the image. Since each wavelet basis function 
V 'ji,h'Juh( x 'y) * s localised at me relevant scale/position, the cor- 
responding wavelet coefficient bj l j 2 -j 1 i 2 measures the amount of 
signal in the image at this scale and position. 

2.1 Orientation sensitivity 

At this point, it is important to note the sensitivity of the orthogonal 
wavelet decomposition to the orientation of the original input map. 
As shown by Daubechies (1992), it is impossible to construct an or- 
thogonal wavelet basis, in which the basis functions are both sym- 
metric (or anti-symmetric) and have compact support. This asym- 
metry of the basis functions is the cause of the orientation sensi- 
tivity. This is most easily appreciated by considering an input map 
consisting of just one of the wavelet basis functions. If this map 
is rotated through 180 degrees (say), then because the basis func- 
tions are asymmetric it is not possible to represent the rotated basis 
functions in terms of just one of the original basis function. In- 
stead, the signal in the rotated map must be represented by several 
wavelet basis functions with different scale and position indices. 
Thus any statistics based on the wavelet coefficients are sensitive 
to the orientation of the original input map. Since the origin of this 
effect is the asymmetry of the one-dimensional wavelet basis func- 
tions, it also occurs for two-dimensional orthogonal wavelet de- 
compositions based on the Mallat algorithm (Mallat 1989), which 
is also commonly called the multiresolution analysis method. In 
order to obtain wavelet statistics that are invariant under 90, 180, 
270 degrees rotations of the input image (and also insensitive to 
cyclic translations of the image by an arbitrary number of pixels in 
each direction), it is necessary to use the a trous wavelet algorithm 
(see e.g. Starck, Murtagh & Bijaoui 1998) with a symmetric filter 
function. The application of this technique to the detection of non- 
Gaussianity in the CMB will be presented in a forthcoming paper. 
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Figure 1. The regions of the two-dimensional wavelet domain for a 32 x 32 
image. The italic numbers show the value of k = 2 ;i + 2' 1 in each region; 
see text for details. 



2.2 Application to COBE data 

In this paper, we will be concerned with Face and Face 5 of 
the COBE QuadCube pixelisation scheme in Galactic coordinates, 
each of which consists of 32 x 32 equal-area pixels (i.e. J\=J2= 5) 
of size (2.8°) 2 . Thus the scale j corresponds to an angular scale of 
2.8 x 2 4 ~- / . Following the discussion by HJL, the structure of the 
corresponding wavelet domain is shown in Fig. [IJ where the pixel 
numbers are plotted on a logarithmic scale. We see that the domain 
is partitioned into separate regions according to the scale indices j\ 
and j2 in the horizontal and vertical directions respectively. Thus 
regions with j\ = 72 contain wavelets basis functions that represent 
the image at the same scale in the two directions, whereas regions 
with j] ^ 72 describe the image on different scales in the two di- 
rections. As discussed in HJL, regions with j\ = or 72 = ac- 
tually contain basis functions that are tensor products of different 
one-dimensional basis and so for the remainder of this paper we 
restrict our attention to regions with 71,72 ^ 1- We also define the 
integer variable k = 2 J> +2 J2 , which serves as a measure of inverse 
scale length, and is constant within each region of the wavelet do- 
main. We note that the value of k is not altered if the values j\ and 
72 are interchanged. In this paper, we also restrict ourselves to the 
Daubechies 4 wavelet basis used by PVF, although analogous anal- 
yses may also be performed for other orthogonal discrete wavelet 
bases, and indeed similar results to those presented in Section 3 are 
obtained. 



2.3 Skewness and kurtosis spectra 

Following HJL, when considering the statistics of the wavelet co- 
efficients bj l j 2 -j l i 2 of an image, it is useful to consider separately 
all those coefficients that share each value of k. For each value of 
k, we then use the corresponding wavelet coefficients to calculate 
estimators of the skewness 5 and (excess) kurtosis K of the parent 
distribution from which the coefficients were drawn. We therefore 
obtain the skewness and kurtosis 'spectra' S(k) and K(k) for the 
image. 

As mentioned in the Introduction, at this point our method 
diverges from that used by PVF in two ways. Firstly, PVF only 
consider regions of the wavelet domain for which j\ — 72 and 
7l 5 72 ^ 1 0- e - k = 4,8,16,32), whereas we consider all regions 
with 71,72 ^ 1- Secondly, we calculate the estimators 5 and K in a 



different way from that adopted in PVF, as follows. At each value 
of k the skewness and (excess) kurtosis of the parent distribution of 
the wavelet coefficients are given by 



„ , 3/2 , 3/2 

S = Pi/p 2 = K 3/V . 

K = J U4//"2-3 = Kt/ K 2 ) 



(3) 
(4) 



where p„ is the nth central moment of the distribution and K„ is the 
nth cumulant (see HJL for a brief discussion). In PVF, the estima- 
tors p„ of the central moments are simply taken to be the central 
moments of the sample of wavelet coefficients. It is easily shown, 
however, that these estimators are biased, so that (p„) 7^ p n , and 
this bias is quite pronounced when the sample size is small (as it 
is in this case). PVF then estimate the skewness and (excess) kur- 
tosis by inserting the biased estimators p„ into (^ and (Q) respec- 
tively. Thus, the corresponding estimators S and K are also signif- 
icantly biased. In this paper, we instead calculate our estimates of 
the skewness and (excess) kurtosis using ^-statistics (see Kenney 
& Keeping 1954; Stuart & Ord 1994; HJL). These provide unbi- 
ased estimates K„ of the cumulants of the parent population from 
which the wavelet coefficients were drawn. These unbiased estima- 
tors of the cumulants are then inserted into ^ and (Q) to obtain the 
estimators 5 and K. 



2.4 Scale-scale correlation spectrum 

In addition to the skewness and kurtosis spectra, we may also mea- 
sure the correlation between the different domains of the wavelet 
transform by defining the estimators of the scale-scale correlation 
as 



C p ■ 

J2 



2/.+A+2£ £ b P 



h 1 ' 



ji,h;[h/*Uh/2] ji+i,h+i-M 



l 2 D h,j 1 ;[l l /2],[l 2 /2] 



Li, £z 2 



(5) 



7i+i,72+i;/i ,h 



In equation (^|), the sums on /; extend from to 2 /l + 1 — 1 (simi- 
larly for I2), p is an even integer and [ ] denotes the integer part. 

Thus C; ; measures the correlation between the wavelet coeffi- 

J1J2 

cients in the domains (71,72) and (71 + 1, 72 + !)■ In PVF, it was 
assumed that j\ = 72, so that the correlation of wavelet coefficients 
were only calculated between adjacent diagonal domains in Fig. ll 
When 71 7^ 72, however, it is convenient to extend the sums in (f) 
to include also the corresponding domains with j\ and 72 inter- 
changed. Thus, in each case, we in fact measure the correlation 
between wavelet coefficients with inverse scalelengths of k and 2k 
respectively (see Fig. [j]). For each possible value of k, we denote 
this correlation by C p (k), thereby producing a scale-scale correla- 
tion spectrum. Following PVF, we restrict our analysis to the case 
where p = 2. 



2.5 The non-Gaussianity test 

The skewness, (excess) kurtosis and scale-scale correlation spectra 
S(k) , K(k) and C 2 (k) of the wavelet coefficients form the basis 
of the non-Gaussianity test. The procedure is as follows. We first 
calculate the S(k), K(k) and C 2 (k) spectra for Face or Face 5 
of the 4-year COBE map. We then generate 5000 realisations of 
an all-sky CMB map drawn from an inflationary/CDM model with 
parameters £l m = 1, £1^ = 0, h = 0.5, n = 1 and 2rms-ps = 18 pK, 
convolved with a 7°-FWHM Gaussian beam. For each realisation, 
we then add random Gaussian pixel noise, where the rms of the 
noise in each pixel is taken from the COBE rms noise map. The 
S(k), K(k) and C 2 (k) spectra are then calculated for Face and 
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Face 5 of each of the 5000 realisations to obtain approximate prob- 
ability distributions for the S{k), K(k) and C 2 (k) statistics when 
the CMB signal is the chosen Gaussian inflationary/CDM model. 
By comparing these probability distributions with the correspond- 
ing spectra for Face and Face 5 of the COBE map, we thus obtain 
(at each £-value) an estimate of the probability that the CMB signal 
in the DMR-DSMB map is drawn from a Gaussian ensemble char- 
acterised by the chosen inflationary/CDM model. For each face, 
however, we obtain skewness and kurtosis statistics at ten different 
^-values, and six different scale-scale correlation statistics. Thus, 
the total number of statistics obtained for each face is 26, and care 
must be taken in assessing the significance of non-Gaussianity de- 
tections at individual ^-values (see below). As discussed in section 



2. 1 , however, the orthogonal wavelet decomposition is sensitive to 
the orientation of the input map. Thus, we repeat the above non- 
Gaussianity test for the case where Face and Face 5 are both ro- 
tated through 180 degrees. 

It is also clear that, to some extent, the results of such an 
analysis will depend on the chosen parameters in the inflation- 
ary/CDM model via the corresponding predicted ensemble-average 
power spectrum Q, from which the 5000 realisations are gener- 
ated. Nevertheless, since at each ^-value the skewness and kurtosis 
statistics contain the variance /./2 of the wavelet coefficients in their 
denominators, and the scale-scale correlation in ^ is similarly 
normalised, we would expect these statistics to be relatively un- 
affected by changing the power spectrum of the inflationary/CDM 
model. As an interesting test, we repeated our entire analysis for 
the case where the 5000 realisations were instead generated using 
the maximum-likelihood Q spectrum calculated from the 4-year 
COBE data by Gorski (1997). As expected, we found that the re- 
sults were virtually identical to those presented in the next Section. 



3 RESULTS 

3.1 The DMR-DSMB map 

In this Section, we present the results of the wavelet non- 
Gaussianity test when applied to Face and Face 5 of the 4-year 
COBE DMR-DSMB map in Galactic coordinates. As mentioned 
in the Introduction, this Galaxy-removed map is constructed by 
first subtracting templates of synchrotron and dust emission and 
then removing the free-free emission (Bennett et al. 1994). We find 
that the results of the non-Gaussianity test are similar for both the 
DMR-DSMB and DMR-DCMB Galaxy-removed maps. 

The resulting S(k), K(k) and C 2 {k) spectra for Face and Face 
5 of the DSMB map are plotted in Fig. [2 In each plot, the crosses 
correspond to the values derived from the DSMB map orientated 
in the same manner as that used by PVF (orientation A), the solid 
squares correspond to the values obtained from the DSMB map af- 
ter rotating it through 180 degrees (orientation B), and the open 
circles denote the mean of corresponding distribution derived from 
the simulated COBE observations of the 5000 realisations of the in- 
flationary/CDM model. The error bars denote the 68, 95 and 99 per 
cent limits of the distributions. These distributions were found to 
be virtually indistinguishable for the two orientations of the COBE 
data. For convenience, the S(k) and K(k) spectra have been nor- 
malised so that the variance of each distribution is equal to unity. 
Thus, for any particular k- value, a estimate of the significance level 
can be read off directly from the scale on the vertical axis. 

As mentioned above, we calculate the S(k) and K(k) spectra 
for all available domains of the wavelet transform, and the C 2 (k) 



spectrum for all pairs of domains whose k- values differ by a factor 
of 2 (with ji , ]2 1 in each case; see Fig. |l|). In contrast, PVF only 
considered domains with j\ = jj and thus only obtained S(k) and 
K(k) values for k = 4, 8, 16, 32, and C 2 (k) values at k = 4, 8, 16. 

We see from Fig. ^ that for orientation A (crosses), all the 
points in the skewness and kurtosis spectra lie comfortably within 
their respective Gaussian probability distributions for both faces. 
In the scale-scale correlation spectrum, however, we confirm PVFs 
finding of a point at k — 4 that lies slightly outside the 99 per 
cent confidence limit. On the other hand, for orientation B (solid 
squares) we obtain two skewness detections someway beyond the 
99 per cent confidence limit. These occur in Face at k = 32 and 
in Face 5 at k = 24. From Fig. [l[ however, we see that these k- 
values correspond to wavelet basis functions on small scales, cor- 
responding to pixel-to-pixel variations in the COBE map. Thus it 
is unlikely that this non-Gaussianity is cosmological in origin; we 
return to this point below. The kurtosis spectrum and scale-scale 
correlation spectra show no strong non-Gaussian outliers for this 
orientation. 

In order to investigate the robustness of the high-A: outliers in 
the S(k) spectra for orientation B, we repeated the analysis for the 
DSMB map with all multipoles above I — 40 removed. A similar 
filtering process was also performed on each of the 5000 CDM re- 
alisations. Since the 7-degree FWHM COBE beam essentially fil- 
ters out all modes beyond I = 40, we would expect these modes 
to contain no contribution from the sky and consist only of instru- 
mental noise or observational artefacts. We also repeated the filter- 
ing process for orientation A. The corresponding S(k), K(k) and 
C 2 {k) spectra for two orientations of Face and Face 5 of the fil- 
tered DSMB map are plotted in Fig. |5| We see immediately that 
the high-/: skewness detections that were present in the unfiltered 
map have now disappeared. This suggests that the non-Gaussianity 
present in the original DSMB map is not cosmological in origin, 
and is most likely an artefact resulting from the algorithm used to 
subtract Galactic emission. From Fig. ^(c), we also note that the 
three points that lay outside the 95 per cent limit in the C 2 (k) spec- 
trum for Face of the original DSMB map in orientation B (see 
Fig.^(c)) have all now been brought well within the Gaussian error 
bars. Thus we find no strong evidence for non-Gaussianity in the 
filtered DSMB map in orientation B. For orientation A, however, 
as we might expect, the level of significance of the C 2 detection at 
k = 4 was only slightly reduced by the filtering process. 



3.2 The 53+90 GHz coadded map 

Since the above analysis suggests some non-Gaussianity on pixel 
scales in the DSMB map, possibly introduced by the Galaxy sub- 
traction algorithm, we repeat the analysis for the inverse noise vari- 
ance weighted average of the 53A, 53B, 90A and 90B COBE DMR 
channels. 

Fig. ^ shows the S(k), K(k) and C 2 (k) spectra for Face and 
Face 5 of the 53 + 90 GHz coadded map in both orientations. For 
orientation A (crosses), none of the skewness, kurtosis or scale- 
scale correlation statistics lies outside the corresponding 99 per cent 
limit. Also, for orientation B (solid squares), we see that, in contrast 
to the DSMB map, no large detections of non-Gaussianity are ob- 
tained at high k in the skewness spectra. Nevertheless, outliers do 
occur at the 99 per cent level for Face in the K(k) spectrum at 
k = 4, and for Face 5 in the C 2 {k) at k = 6 and k = 12. Indeed, 
the last of these lies someway outside the 99 per cent confidence 
limit. However, this statistic measures the correlation between the 
wavelet coefficients in the domains with k = 12 and k = 24, and is 
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Figure 2. The S(yfc), K(k) and C 2 (<t) spectra for Face [plots (a), (b) and (c)] and Face 5 [plots (d), (e) and (f)] of the DMR-DSMB map. The crosses correspond 
to the orientation used by PVF (orientation A), whereas the solid squares correspond to the map rotated through 180 degrees (orientation B). The open circles 
show the mean value for each statistic obtained from 5000 simulated COBE observations of CDM realisations, and the error bars denote the 68, 95 and 99 per 
cent limits of the corresponding distribution. For convenience, the S(k) and K{k) spectra have been normalised at each value of k such that the variance of the 
distribution obtained from the 5000 CDM realisations is equal to unity. 
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Figure 3. As in Fig. |, but for the DMR-DSMB map with all mutlipoles above £ = 40 removed. 



therefore influenced primarily by features in the map on the scale 
of one or two pixels in size. 

We once again tested the robustness of these putative detec- 
tions of non-Gaussianity by repeating the analysis after removing 
all multipoles above £ = 40 from the COBE map and the CDM re- 
alisations. The resulting spectra are shown in Fig. [5]. From Fig. ^|(f), 
we see that the large outlier in C 2 (12) that was obtained for the un- 
filtered map in orientation B has now reduced to well within the 
Gaussian error bars. This suggests that the noise in Face 5 of the 



coadded map may contain some non-Gaussian component. Never- 
theless, the two outliers at the 99 per cent limit in K(4) for Face 
and C 2 (6) for Face 5 in orientation B remain unaffected by the 
filtering process, and thus might be interpreted as robust signatures 
of non-Gaussianity on large scales. 

It is, however, important to remember that, although the sig- 
nificance level is above the 99 per cent level for these individual 
statistics, we must take into account the fact that no outliers are 
found in the large number of other statistics we have calculated; 
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Figure 4. As in Fig. Q but for the 53+90 GHz coadded map. 
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Figure 5. As in Fig. ^, but for the 53+90 GHz coadded map with all multipoles above £ = 40 removed. 



this is discussed below. It should also be bourne in mind that no 
Galaxy subtraction has been performed on the 53+90 GHz coad- 
ded map. Although, our analysis is restricted to Face and Face 
5 of the COBE QuadCube, which lie outside the standard Galac- 
tic cut, it is possible that these faces may be contaminated to some 
extent by high-latitude Galactic emission. 



4 DISCUSSION AND CONCLUSIONS 

We have presented an orthogonal wavelet analysis of the 4-year 
COBE data, in order to search for evidence of large-scale non- 



Gaussianity in the CMB. In particular, we identify an orientation 
sensitivity associated with this method, which must be borne in 
mind when assessing its results. 

We find that several statistics in the S(k), K(k) and C 2 {k) spec- 
tra for the COBE DSMB and 53+90 GHz coadded maps (in orienta- 
tions A and B) lay outside the 99 per cent limit of the corresponding 
probability distributions derived from 5000 simulated COBE obser- 
vations of CDM realisations. However, only one such outlier in the 
DSMB map and two outliers in the 53+90 GHz coadded map were 
found to be robust to the removal of all multipoles above I = 40 
in the COBE map and CDM realisations. In the DSMB map, this 
occurs in C 2 (4) for Face in orientation A, and in the 53+90 GHz 
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coadded COBE map the outliers are in K(4) for Face and C 2 (6) 
for Face 5, both in orientation B. 

We must, however, take care in assessing the significance of 
these outliers. For each face and orientation we calculate 26 differ- 
ent statistics. Thus for each data set (either DSMB or 53+90 GHz 
coadded), the total number of statistics is 2 x 2 x 26 = 104, and we 
must take proper account of the fact that a large number of these 
show no evidence of non-Gaussianity (see, for example, Bromley 
& Tegmark 1999). Since the statistics presented here are not inde- 
pendent of one another and generally do not possess Gaussian one- 
point functions, the only way of obtaining a meaningful estimate 
of the significance of our results is by Monte-Carlo simulation. In- 
deed, in their bi-spectrum analysis of the 4-year COBE data, Fer- 
reira et al. (1998) used Monte-Carlo simulations and a generalised 
% 2 -statistic to assess their results. In our case, we adopt a slightly 
different approach and simply use the 5000 CDM realisations to 
estimate the probability of obtaining a given number of robust out- 
liers at > 99 percent level in any of our 104 statistics, even when 
the underlying CMB signal is Gaussian. For the DSMB data, we 
obtained one robust outlier, and the corresponding probability of 
this occuring by chance was found to be 0.59. For the 53+90 GHz 
coadded data, two outliers were obtained, and the corresponding 
probability is 0.28. Therefore planar orthogonal wavelet analysis 
of the 4-year COBE data can only rule out Gaussianity at the 41 
per cent level in the DSMB data and at the 72 per cent level in the 
53+90 GHz coadded data. Thus, we conclude that this method does 
not provide strong evidence for non-Gaussianity in the CMB. 



Pando J., Valls-Gabaud D., Fang L.-Z., 1998, Phys. Rev. Lett., 81, 4568 
(PVF) 

Press W.H., Teukolsky S.A., Vettering W.T., Flannery B.P., 1994, Numerical 

Recipes. Cambridge University Press, Cambridge 
Schmalzing J., Gorski K.M., 1998, MNRAS, 297, 355 
Starck J.-L., Murtagh F., Bijaoui A., 1998, Image Processing and Data 

Analysis. Cambridge University Press, Cambridge 
Stuart A., Ord K.J., 1994, Kendall's Advanced Theory of Statistics (Volume 

I). Edward Arnold, London 
Turok N., 1996, ApJ, 473, L5 

This paper has been produced using the Royal Astronomical Soci- 
ety/Blackwell Science DTpX style file. 



ACKNOWLEDGEMENTS 

The authors thank David Valls-Gabaud for his work in indepen- 
dently verifying the numerical results presented here. PM acknowl- 
edges financial support from the Cambridge Commonwealth Trust. 
MPH thanks the PPARC for financial support in the form of an 
Advanced Fellowship. 



REFERENCES 

Banday A.J., Gorski K.M., Bennett C.L., Hinshaw G, Kogut A., 
Lineweaver C, Smoot G.F., Tenorio L., 1997, ApJ, 475, 393 

Banday A.J., Zaroubi S.. Gorski K.M.. 2000. ApJ. in press (a stro- 
ph/9908070) 

Bennet C.L. et al., 1992, ApJ, 396, L7 

Bennet C.L. et al., 1994, ApJ, 436, 423 

Bennet C.L. et al., 1996, ApJ, 464, 1 

Bouchet F.R., Bennett D.P, Stebbins A., 1988, Nat, 335, 410 
Bromley C.L., Tegmark M., 1999, ApJ, 524, L79 
Colley W.N., Gott III J.R., Park C. 1996, MNRAS, 281, L82 
Daubechies I., 1992, Ten Lectures on Wavelets. S.I. A.M., Philadelphia 
Ferreira P., Magueijo J., Gorski K.M., 1998, ApJ, 503, LI 
Gorski K.M., 1997, in Bouchet F.R., Gispert R., Guideroni B., Tran Thanh 
Van J.,eds, Proc. 16th Moriond Astrophysics Meeting, Microwave 
Anisotropics. Editions Frontiere, Gif-sur-Yvette, p. 77 
Heavens A.F, 1998, MNRAS, 299, 805 

Hobson M.P., Jones A.W., Lasenby A.N., 1999, MNRAS, 309, 125 (HJL) 
Kenney J.F., Keeping E.S., 1954, Mathematics of Statistics (Part I). Van 

Nostrand, New York 
Kogut A., Banday A.J., Bennett C.L., Gorski K.M., Hinshaw G., Smoot 

G.F., Wright E.L., 1996, ApJ. 4 64.1.29 

Magueijo J., 1999, ApJ, submitted ( ^tstro-ph/991 133^ ) 

Mallat S.G., 1989, IEEE Trans. Acoust. Speech Sig. Proc, 37, 2091 

Novikov D., Feldman H.A.. Shandarin S.F., 1999, Int. J. Mod. Phys. D, 8, 

291 



© 1998 RAS, MNRAS 000, 



